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PACS numbers: 

We present scanning SQUID microscopy data on the 
superconductors Sr2Ru04 (Tc = 1.5 K) and PrOs4Sbi2 
(Tc = 1.8 K). In both of these materials, superconductivity- 
related time-reversal symmetry-breaking fields have been 
observed by muon spin rotation; our aim was to visu- 
alize the structure of these fields. However in neither 
Sr2Ru04 nor PrOs4Sbi2 do we observe spontaneous 
superconductivity-related magnetization. In Sr2Ru04, 
many experimental results have been interpreted on the 
basis of a ± ipy superconducting order parameter. This 
order parameter is expected to give spontaneous magnetic 
induction at sample edges and order parameter domain 
walls. Supposing large domains, our data restrict domain 
wall and edge fields to no more than ~0.1% and ~0.2% of 
the expected magnitude, respectively. Alternatively, if the 
magnetization is of the expected order, the typical domain 
size is limited to ~30 nm for random domains, or '^500 nm 
for periodic domains. 

Introduction 

Sr2Ru04 is a highly two-dimensional, layered per- 
ovskite superconductor with Tc ~ 1.5 K in the clean 
limit. Its in-plane coherence length is iah(S^) = 66 nm [1], 
and magnetic penetration depth \ah = 160-190 nm [2]- 
[4j On the Fermi sheet thought to be dominant for su- 
perconductivity (the 7-sheet) the electron mass is 16me 
(5.5 times the band mass) [l]. PrOs4Sbi2 is a heavy 
fermion superconductor of cubic symmetry, with m* in 
the range of tens of rrie [5l . Measurements of specific 
heat show an unusual double superconducting transition, 
with Tci ^ 1.86 K and Tc2 ^ 1.70 K [7 ,9j. The coherence 
length of PrOs4Sbi2 is ^ = 12 nm, while muon spin rota- 
tion (/xSR) measurements yield A(0) ^ 350 nm p!Q[ ITT]. 

Sr2Ru04 and PrOs4Sbi2 both have unusual supercon- 
ducting states. Strong evidence for either line nodes 
or deep gap minima, or both, in Sr2Ru04 comes from, 
e.g.^ RF [12 , microwave [2 and specific heat measure- 



ments [T3l[T4]. In PrOs4Sbi2, RF and thermal conduc- 
tivity measurements suggest point nodes [15 , while /iSR 
and Sb NQR show fully-gapped superconductivity (at 
ambient pressure) JOl Ull USj. Nodes may exist on a 
small-gap band suppressed by modest applied fields [11] . 
Additional thermal conductivity measurements confirm 
multi-band superconductivity, but indicate no nodes on 
either gap p/7|. In Sr2Ru04, specific heat measured 
against applied field confirms multi-band superconduc- 
tivity [T3], and that if not nodes there are at least deep 
minima on the primary band. 

Triplet pairing has been shown in Sr2Ru04 with high 
certainty by measurement of the ^^O Knight shift |T8j. 
Measurement of the muon Knight shift in PrOs4Sbi2 also 
indicates triplet pairing, but with less certainty: owing to 
low-lying crystal electric field states, the expected Knight 
shift for singlet pairing is less clear [19]. 

The absence of a Hebel-Slichter peak in Ru NMR mea- 
surements on Sr2Ru04 [20], and in Sb NQR measure- 
ments on PrOs4Sbi2 p!6|[2T^ indicate unconventional su- 
perconductivity in both materials. In Sr2Ru04 it is con- 
firmed by demonstration that Tc ^ as the mean free 
path shrinks to ^iah f22]. In PrOs4Sbi2, Td shows a 
modest, and Tc2 possibly a more pronounced, sensitiv- 
ity to sample quality [7 . Odd-parity orbital symme- 
try in Sr2Ru04 has been shown by fabrication of a tt- 
SQUID [23], and superconductivity-related time-reversal 
symmetry-breaking (TRSB) by measurement of the Kerr 
effect [24]. A two-component order parameter is indi- 
cated by Josephson interferometry [25], hysteretic trans- 
port in microstructures [26], and a jump in the transverse 
sound velocity at Tc [27]. These results have all been in- 
terpreted in terms of a chiral Px i Wy orbital order pa- 
rameter. The order parameter of PrOs4Sbi2 remains an 
open question; possibilities for both singlet and triplet 
pairing are listed in Ref. [28] . 

Along with the shared features described above, 
Sr2Ru04 and PrOs4Sbi2 are also the only two materi- 
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als where observation of spontaneous, superconductivity- 
related TRSB fields by /iSR is well-established. In 
Sr2Ru04, an average internal induction far below Tc 
of ~0.5 G is found, with the rapid initial decay of 
muon polarization indicating a peak induction of at least 
~5 G [4] [29l, indicating a dilute density of sources. In 
PrOs4Sbi2, the average internal induction is at least 
twice as large, ~1.5 G [30 , but its distribution is not as 
peaked as in Sr2Ru04, suggesting a higher source den- 
sity. The TRSB appears to onset at the upper (^1.85 K) 
Tc. 

For px i Wy order in Sr2Ru04 a magnetization along 
the z axis (crystalline c axis) is expected: the orbital 
angular momentum of the condensate would give an un- 
cancelled current at domain edges (meaning domain walls 
and sample edges). Inward from these edges, Meissner 
screening would in turn result in counterfiowing screen- 
ing currents. If each pair in the 7-sheet condensate 
is assigned angular momentum fi^ an edge current of 
kphe/87rm* =2.6 fiA per layer results jSj, for a field 
discontinuity of 50 G. Matsumoto and Sigrist (MS) have 
solved the Bogoliubov-de Gennes equations in a quasi- 
classical approximation for an ideal Px ± ipy supercon- 
ductor (without secondary bands or gap minima/nodes, 
and with specular scattering at edges) and obtain edge 
and domain wall inductions peaking at ~10 and ~20 G, 
respectively [32 . (We scale their unitless results by 
^ = 66 nm and = 165 nm; they assume = 2.5.) 

There is confiicting experimental guidance on domain 
size: Kerr rotation indicates domains at least a few 
times larger than the beam size, or ^100 /im [24 , while 
Josephson interferometry suggests domains at the edges 
~1 jam across [25]. The tt- SQUID required phase coher- 
ence across the 0.6 mm width of the Sr2Ru04 crystal, 
indicating large domains [23|. If the MS result is ap- 
proximately correct {i.e., a ~10 G field across a width 
~2Aa6 at domain walls), the 0.5 G average interal field 
observed by /iSR suggests ~10 jam domains. Domain 
size might be affected by the cooling rate through Tc, 
which was ~1 K/hr for the Kerr, Josephson interferom- 
etry, and TT-SQUID experiments, and faster for the /iSR 
measurement. Domain size is discussed in more detail in 
Ref. [33]. 

To date, edge and domain wall fields have not been ob- 
served by scanning magnetic probes in Sr2Ru04. Scan- 
ning Hall probe measurements by Bjornsson et al con- 
strain edge and domain wall currents to be less than ~3 
and 8% of the MS results, repectively [34l [35] [The scan 
area was (70 /im)^, so domain walls may have been ab- 
sent.] Kirtley et al, in SQUID scans spanning a ~1 mm- 
wide sample, improve the limit on both edge and domain 
wall currents to about 1% of the expectation, for domains 
larger than ~8 /im [35^ . 

Analysis by Ashby and Kallin show that nonspecular or 
pair-breaking edge scattering could reduce the expected 
edge magnetization, but not by the orders of magnitude 
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FIG. 1: Left: photograph of the Sr2Ru04 sample scanned 
in this work. The surface is a cleaved ab surface. The left 
edges are polished and the right is a growth edge. Cleave 
terraces visible in the magnetic scans are circled. Middle: 
PrOs4Sbi2. The right edge is a growth edge. Right: front end 
of the SQUID used in this work. The pick-up coil (of diameter 
3.2 /im) and contact tip, the point where the SQUID contacts 
the sample, are indicated. 

required for consistency with experiment [36]. Selection 
of Ginzberg-Landau parameters nearer the edge of sta- 
bility for px + ipy order reduces the edge currents, but 
also spreads them out over a larger range, so the total 
fiux, and the observed signal in scanning probes microns 
above the surface, would not be greatly reduced. 

The work here further tightens the limits on chiral 
currents and domain structure in Sr2Ru04. We also 
discuss the possibility of periodic domains, and show the 
first magnetic scans of PrOs4Sbi2. 

Magnetic scans of Sr2Ru04 

The scanning SQUID used here is a niobium-based de- 
vice, described in Refs. [3'^ and [38]. Flux is coupled 
into the SQUID through a 3.2 /im-diameter pick-up coil 
(Fig. [1]). The leads to the pick-up coil are shielded to 
minimize fiux coupling into the space between the leads. 
SQUIDs are fiux-sensitive devices, so the units on the 
data shown in this work are units of fiux; 1 <I>o = hc/2e 
in a 3.2 jiia loop corresponds to an average induction of 
^2.5 G. The SQUID was generally scanned on a plane 
~1 jiia out of contact with the sample, to reduce noise 
and spurious features from surface roughness. 

The Sr2Ru04 crystal (photograph: Fig. [T]) was grown 
in a fioating zone furnace, and was not annealed after 
growth [39 . All scans were done at T = 0.4 K, with fairly 
rapid cooling through Tc, ~1 K/min. For the mosaic of 
scans shown in Fig. [2j the crystal was cooled in Earth's 
~l/2 Oe field with the out-of-plane component cancelled 
by ~50% by an applied field, however the applied field 
was subsequently turned off and vortices reintroduced 
into the sample, to approximately match the ambient 
field, by electromagnetic noise from the positioners. The 
weak tails that extend leftwards from each vortex, clearly 
visible in the close-up of a single vortex, are artifacts of 
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FIG. 2: Three overlapping SQUID magnetometry scans of the 
Sr2Ru04 crystal pictured in Fig.[l] in Earth's ^1/2 Oe field. 
Single vortices are clearly visible. Across the entire scan range 
vortices are lined up along upper-left-to- lower-right stripes. 
A close-up of a single vortex is shown at the upper right; the 
weak tail extending leftwards from the vortex is an artifact of 
the SQUID imaging kernel. 



the imaging kernel: the shielding of the pick-up coil leads 
is not perfect. The full extent of these tails can be seen 
in Fig. [ojc): ^70 /im, about the distance between the 
pick-up coil and the point where the shielding becomes a 
fully-formed coaxial cable. 

Although not the main point of this paper, the mixed- 
state scans contain three features worth noting. (1) The 
local vortex distribution is uneven. Even at small ap- 
plied fields, where direct vortex- vortex interaction is neg- 
ligible, minimization of global field energy encourages a 
homogeneous vortex density — compare with the much 
more homogeneous vortex distribution in clean areas of 
PrOs4Sbi2, at an applied field of 760 mOe [Fig. [6ja)]. 
Local vortex coalescence has been previously reported 
in Sr2Ru04 [Ml SQI E], and also in MgB2, where it 
was explained as originating from the two-component 
order parameter, one in the type-II regime, and the 
other type-I [42 . As Sr2Ru04 is only weakly type- 
II {Kab = Kb{0)/U{0) ^ 2.6; K > 1/V2 indicates 
type-II) and has a two-component order parameter, the 
cause of local clustering might be similar. (2) The large- 
scale distribution is also uneven: towards the lower right 
and upper left, vortex- free areas up to ~30 fim across 
are adjacent to similar-sized regions of high vortex den- 
sity. These regions may indicate spatially-varying sam- 
ple quality. (3) The distribution is anisotropic: over the 
entire ~300 /im-wide scan area, vortices line up along 



upper-left-to- lower-right stripes. A striped vortex distri- 
bution was also reported in [34 . In highly anisotropic 
superconductors c-axis vortices can form chains along in- 
plane vortices, a phenomenon that has been observed 
in Sr2Ru04 [43^. However an in-plane field of ~10 Oe 
appears to be required to form and orient well-defined 
chains [40], whereas the field for Fig. [2] was no more than 
Earth's -1/2 Oe, and -50 mOe in Ref. [34 . 

Fig. [3] shows this paper's main results on Sr2Ru04: 
magnetic scans in the near-absence of vortices. The 
scans overlap, however each is a separate thermal cycle 
(to above Tc). Four prominent features appear in the 
scans: (1) In panel (a), there are two vortices, one at 
the lower left, and the other beyond the right edge of 
the scan. (The extended dipole-like feature is part of its 
tail.) (2) The step changes in the signal across cleave 
terraces and the sample edge result from SQUID-sample 
interaction, not static magnetization: the SQUID is 
scanned on a plane —1 jam above the sample surface, so 
when the pick-up coil passes over topographic features 
its proximity to the sample changes. In operation, the 
SQUID is biased with a DC voltage and the current, 
which varies with the SQUID's critical current, is 
measured. Changing the proximity between the pick-up 
coil and a metallic surface affects the inductance of that 
arm of the SQUID, which in turn affects the critical 
current (as does varying flux in the SQUID, the desired 
signal) [44 . (3) An edge shadow, the blurred dark line 
in panels (b) and (c), is also an artifact. It appears to 
correspond to a step change in magnetic field across 
the sample edge, with the —70 /im distance between 
the shadow and the edge being set by the far extent 
of the SQUID imaging kernel, i.e. the tails discussed 
above. (4) There are numerous magnetic dipoles, most 
very weak: compare their signals with that from the 
vortex. The source of the dipoles was not investigated, 
however typical dipoles in the scans, with peak signals 
of 0.1-5 m<l>o, could be accounted for by inclusions 
of SrRuOa, a ferromagnet with 1.6-2 /i^/Ru [45] . 
—20-80 nm on a side. (Tails on some of the dipoles, like 
the vortex tails, are imaging kernel artifacts.) 

The absence of observed spontaneous mag- 
netization in Sr2Ru04 

No edge or domain wall magnetization, expected for 
Px ± ipy order, is apparent in the Sr2Ru04 scans. To 
compare this null result with theoretical expectation, it 
is necessary to account for the finite sensor resolution 
and scan height. We model the pick-up coil as a 3.2 /im- 
diameter wire loop parallel to the surface, with perfect 
coupling of Bz inside and zero outside. Empirical scan 
heights consistent with this model are obtained by study- 
ing the vortex and dipoles in Fig. [3j For scan heights 
much larger than Aa6, the field distribution of a vortex 
approaches that of a monopole placed a depth Xab be- 
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FIG. 3: Overlapping scans of Sr2Ru04 at 0.4 K; each is a sepa- 
rate thermal cycle. In none of these three scans does the spon- 
taneous edge and domain wall magnetization expected from 
chiral superconductivity appear, (a) A scan over the cleave 
terraces, of total height ^700 nm, indicated in Fig.[T]^a), and 
(lower panel) a section along the dashed line. Note the bro- 
ken color scale: the peak vortex signal far exceeds all other 
features. The terrace signal is an artifact of SQUID-sample 
interaction. The vortex tail, from a vortex ^50 fim beyond 
the right edge of the scan, is an imaging kernel artifact, (b) 
An area with fewer terraces. The "edge shadow" is an imag- 
ing kernel artifact. Other features include dipoles (magnetic 
inclusions) and some terraces, (c) A scan over the sample 
edge, and (lower panel) a section along the dashed line. The 
rms signal in the solid and dotted boxes, after local plane 
subtraction, are 0.06 m$o and 0.05 m$o, respectively. 



neath the surface. The vortex in Fig. (sja) has a FWHM 
(along y) of 3.8 /im and a peak signal of 0.25<l>o. In 
the monopole model, these values correspond to heights 
above the monopole of 1.6 and 1.9 /im, respectively, 
or, subtracting A, ^1.6 /im above the sample surface. 
Heights above the dipoles are obtained by fitting the 2D 
scan data to the simulated response from a point-like 
dipole, with dipole strength, dipole orientation and scan 
height as free parameters. The very prominent dipole 
near the bottom of Fig. [sj^b) gives a height of 1.7 /im. 
Fits to two other dipoles in Fig. |3jb) give heights around 
1.5 /im, and to two dipoles in Fig. jsja), a few 0.1 /im 
higher. In the Sr2Ru04 scans, the scan height above the 
surface is z 1.5 /im. 

The MS calculation gives the magnetic induction deep 
beneath the upper surface of the sample. We extend 
their results to the space above the sample following the 
procedure in Refs. [35 and [46 : If the source inductions 
within the superconductor are along then 
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where K = ^/W^^[^(JJ\ab^ ^ ^^(k) is the Fourier trans- 
form of ^^(x), and z is the height above the sample 
surface. The K/ (k + K) prefactor arises from Meissner 
screening: dX z ^ \^ its effect is that the field distribu- 
tion at height z is what it would be at 2: + A without the 
prefactor. We obtain the expected signal at edges and 
domain walls by applying eq. 1, with A^b = 190 nm, to 
the MS results, then averaging over a 3.2 /im-diameter 
circle. (The accounting of Meissner screening in eq. 1 
is correct over featureless surfaces. However we apply 
eq. 1 to edges, too: the error on z exceeds A, so this is 
a secondary error that does not affect our conclusions.) 
Comparisons of expected and observed signals are shown 
in Fig. [4j The expected signal is shown for z = 1.5 /im, 
however the choice of z is not critical: a doubling of z 
reduces expected edge signals by ~30% and domain wall 
signals by ~60%. 

Fig. [sja) and (b) show simulated scans, taking z = 
2 /im (the upper bound of the actual scan height), for par- 
ticular configurations of domains averaging 15 and 2 /im 
across, respectively. The domains are quasiperiodic for 
computational convenience; their rms area variation is 
38% of the mean. In panel (d) the simulated scans mul- 
tiplied by a scale factor are added to the observed signal. 
For domains 2 /im and larger, domain wall fields would 
have been visible at '^0.1% of the MS result. To look at 
edge fields separately, in (e) the edge fields alone of the 
simulation are added to the data. Weak edge fields ap- 
pear as modulations of the shape of the step, and would 
have been visible at ~0.2% of the expectation for 4 /im 
or larger domains, or ~1% for 2 /im domains. 

We cannot definitively exclude the possibility that the 
entire sample is a single domain, and the edge signal is 
merged with the 6 m<l>o step at the sample edge (due 
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FIG. 4: (a-c) Expected signals compared with observed data 
on Sr2Ru04 for (a) a sample edge, (b) an isolated domain 
wall, and (c), domains periodic along one direction in the ab- 
plane. The expected signals were obtained by extension of the 
MS results to a scan height of 1.5 /xm, using eq. 1, followed 
by averaging over a 3.2 /xm-diameter pick-up coil. The data 
in (a) and (b) is along the dashed line in Fig. [sjc), and the 
noise level indicated in (c) is the pixel-to-pixel noise in the 
boxes in Fig.jsj^c). (d) Hypothetical c-axis domain structures, 
with the magnetization direction and equivalent edge currents 
indicated. Taking |M| = 10 G, for the cases at left and right 
our scans indicate L < 20 and 400 nm, respectively. 



to SQUID-sample interaction). In this case the limit on 
the edge signal is ~6 m<l>o, or ~3% of the expectation. 
However, we note again that the 0.5 G induction observed 
by /iSR is a volume average: it is highly unlikely that its 
source is entirely excluded from our ~(150 /im)^ scan 
area. 

How large could domains be if they were completely 
random? For random domains magnetized along z, the 
expected field at scan height z is given by eq. 15 of 
ref. i46j: 
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where M is the typical magnetization and V typical do- 
main volume. Bz would vary on a length scale ~ z] here 
z is comparable to the pick-up coil radius, so the signal 
is ~ Bzz'^. The observed signal over a featureless area 
of the sample is 0.06 m^o rms [see Fig. [sj^c)], and over 
a section of vacuum, 0.05 m<l>o rms. Taking z = 2 /im, 
obtaining a signal less than 0.02 m<l>o with M = 10 [1] G 
requires V < (6 nm)^ [< (30 nm)^]. 



(a) Simulation: 
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FIG. 5: (a,b) Simulated scans of a quasiperiodic domain 
structure in Sr2Ru04. The MS results are extended to z = 
2.0 /im using eq. 1, then averaged over a 3.2 /xm-diameter pick- 
up coil, (c) Observed signal, from Fig. [sjc). (d,e) Observed 
signal plus a scale factor times the simulation. In (e), domain 
wall contributions are removed. 



Larger domains are not excluded if they are periodic. 
Fig. [4]^c) shows the expected signal from domains pe- 
riodic along one direction in the ah plane, taking the 
MS result for the domain wall magnetization. Domains 
perfectly periodic along a could have been as wide as 
~0.5 /im without being detected. 

Domains could also be periodic along c. To estimate a 
limit on c-axis domain widths, we take M = (0, 0, ±M) 
everywhere, with M constant; i.e. the magnetization 
is equivalent to thin sheet currents at the sample edge. 
Meissner screening does not strongly affect the field dis- 
tribution for domain thicknesses L <C A, and so is ne- 
glected. Two possibilities for c-axis periodic domains are 
illustrated in Fig. |4jd). If the domains are all of equal 
thickness L, then the upper half of the top domain gives a 
far-field \B\ linear in L. (All remaining edge current is in- 
corporated into alternating dipoles.) If \M\ is 10 G, then 
for a 3.2 /im-diameter pick-up coil at heights z = 1, 1.5 
and 2 /im, L = 0.02 /im gives peak edge signals of 2.5, 1.8 
and 1.5 m<l>o, respectively. Alternatively, if the thickness 
of the top domain is halved the long-range field scales as 
L^. For L = 0.40 /im and z = 1, 1.5 and 2 /im, the peak 
edge signal is expected to be ^3.8, 2.3 and 1.5 m<l>o (in 
the absence of Meissner screening, which would reduce 
the signal somewhat further). 

A brief summary on the limits on chiral currents and 
domains sizes in Sr2Ru04: 

1. For domains of Px ± ipy order larger than 2 /im 
across, domain wall currents are at most ~ 0.1% of 
the expectation (as calculated by MS). Edge cur- 
rents are at most ~0.2% of expectation for 4 jam 
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and larger domains, and ~1% for 2 /im domains. 

2. If domain wall currents are of the expected mag- 
nitude and domains are perfectly periodic stripes 
in the ab plane, the domain width can be up to 
~0.5 /im (periodicity 1 /im). 

3. If |M| 10 G [1 G] and the domains are random, 
their volume is at most ~(6 nm)^ [~(30 nm)^]. 

4. If |M| 10 G and M is periodic along c, the do- 
main width can be up to ~ 20 nm if all domains are 
of the same width, or ~ 400 nm if the top domain is 
half the width of the others, to give an edge signal 
< 1 m<l>o. 

If a reason were found for chiral edge and domain wall 
currents to be vastly less than the MS result, then a 
separate explanation for the source of the /iSR signal 
would be required. 
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Magnetic scans of PrOs4Sb 
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As with Sr2Ru04, in anticipation of possible edge cur- 
rents an area of the PrOs4Sbi2 extending over an edge 
was selected for scanning. Scans of the mixed state under 
H < 1 Oe are shown in Fig. [6j [The edge is of a deep 
trench rather than the sample edge; a spread-out vor- 
tex, far beneath the SQUID, appears at the lower right 
of Fig. |6|c).] Also as with Sr2Ru04, the cooling rate 
through Tc was ~ — 1 K/min. 

The vortices clustered strongly at extended sample de- 
fects, the more prominent of which are visible in the pho- 
tograph. Where the sample is relatively free of defects 
[the lower right portion of Fig. [6]^ a)], the vortex distribu- 
tion is homogeneous, in contrast to Sr2Ru04. 

A few opposite-sign vortices appear in panel (b), where 
the cooling field was 160 mOe, and in the same area of a 
300 mOe scan. They do not appear in panel (c), where 
the cooling field was near zero, and so are not due to 
a magnetic inclusion. They may be a consequence of 
trapping of positive-sign vortices combined with overall 
fiux expulsion: the vortex density in the 760 mOe scan 
approximately matches the cooling field, but a significant 
portion of the fiux was expelled in the 160 mOe scan 
(even after accounting for the nonzero background field). 

Like the cleave terraces in Sr2Ru04, surface defects 
appear clearly in the nearly vortex-free scan [Fig. [6]^c)], 
as a results of SQUID-sample interaction. This is proved 
by reversing the SQUID bias: the surface defect signal 
changes sign whereas the edge shadow, resulting from 
static magnetic fields at the edge, does not (Fig. [t]). 

The scan height over the PrOs4Sbi2 sample can be 
determined by studying a few vortices, as in Sr2Ru04. 
The isolated vortex in Fig. [6jc) has a FWHM of 3.8 /im 
and a peak signal of 0.30<l>0 7 giving scan heights (above 
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FIG. 6: (a-c) Magnetic scans of the PrOs4Sbi2 sample shown 
in Fig.[l] All scans are at 0.4 K, under different cooling fields, 
Ha- Vortices cluster at extended defects. The indicated vor- 
tices in (b) are used for scan height determination. In (c), 
Ha = 19 mOe cancels the background z-axis field. Four vor- 
tices remain towards the left side of the scan, and a fifth in a 
deep trench on the right side. Note the broken colour scale: 
the vortex signal greatly exceeds other features. Away from 
the vortices, surface defects appear as an artifact of SQUID- 
sample interaction. A shadow ^70 /im left of the edge is an 
artifact of the imaging kernel and the edge. No magnetic fea- 
tures appear that could explain the ^1.5 G internal induction 
observed by /xSR. After local plane subtraction, the rms sig- 
nal in the solid box is 11 m$o, and in the dotted box, 10 m$o- 
(d) Sections along the solid and dashed lines in (c). 
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FIG. 7: Scans of the same area of PrOs4Sbi2, with the 
SQUID bias current reversed between the two scans. Fea- 
tures resulting from SQUID-sample interaction change sign, 
while features resulting from static magnetic inductions do 
not. 



the monopole) in the monopole model of 1.6 /im and 
1.5 /im, respectively, or, subtracting A, z ^ 1.3 /im. The 
indicated vortices in panel (b) yield z ^ 0.9 and 1.0 /im, 
or ^0.5 /im lower than in the Sr2Ru04 scans. 

At a similar level to Sr2Ru04, no superconductivity- 
related TRSB fields are visible in the PrOs4Sbi2 scans. 
Applying eq. 2, for the case of random magnetic do- 
mains, with M = 1 G, A = 0.35 /im and z = 1 /im, a 
signal of less than 0.1 m^o implies domains of volume 
less than (30 nm)^. The lower scan height allows a 
somewhat tighter limit on periodic domains than in 
Sr2Ru04: if the domain magnetization is comparable 
to the MS result, the upper limit on domain width is 
~0.4 /im. 

Discussion and Conclusions 

Ref. [36 discusses the possibility that domain wall 
fields in line with expectation account for the /iSR data, 
while edge fields are sharply reduced. Two mechanisms 
to achieve this are discussed: pair-breaking edge scatter- 
ing combined with a balancing of Ginzberg-Landau pa- 
rameters, and a competing edge order. If the induction 
observed by /iSR is from domain walls, however, then 
as discussed above it is unlikely that observable domain 
walls would have been absent from the scans here and 
in Ref. [35]. Mechanisms that might eliminate edge and 
domain wall fields, while maintaining px ± ipy order, in- 
clude a balancing of GL parameters such that the Px and 
Py components are everywhere nearly equal [47], or the 
effects of multiple bands, or nodes / deep gap minima, 
or strong spin-orbit coupling [48]. However the origin 
of the TRSB fields observed in /iSR would be left un- 
explained. Ref. [30] discusses, in regards to PrOs4Sbi2, 
the possibility of a finite hyperfine field induced at the 
/i+ sites, which would require a nonunitary order param- 
eter. Nonunitary pairing under H = has not been 
confirmed in any material, however, and is considered 
unlikely for Sr2Ru04 pi. Magnetization through pair- 



breaking at impurity sites is not likely to be the source 
of the /iSR signal: observable magnetization did not ap- 
pear at defect sites, of which there must have been very 
many, across a ~100 jam range, in the scans here. Also, 
/iSR data on a lower- Tc sample, with more defects, indi- 
cate weaker, not stronger, TRSB fields [4]. 

As described above, random static domains of magne- 
tization consistent with the /iSR data would need to be 
on a scale ^ or smaller to have evaded detection. It 
seems unlikely that domains of an orbital order param- 
eter could be so small. Also, in Sr2Ru04, the field dis- 
tribution observed by /iSR is more consistent with dilute 
sources than dense random magnetization. 

The limits on periodic domains are less severe. If edge 
currents are of the predicted order (accounting for the 
/iSR data), then the most plausible scenario described 
here for eliminating long-range edge fields may be do- 
main periodicity along the c-axis. If the periodicity is 
good then the domains could be up to ~400 nm ^lOO^c 
in width, implying an energy cost below 1% of the con- 
densation energy. 

Closely-spaced c-axis domain walls extending through- 
out the sample are not ruled out by experiments to date: 
if the structure of a domain wall is such that one of kx 
or ky remains finite across the wall, then by slow cool- 
ing it may be possible to obtain one of kx or ky finite 
throughout the crystal — domain formation may be con- 
trolled by dynamics near Tc, where ^ is longer and do- 
main walls may interact more strongly — which would 
permit fabrication of tt-SQUIDs [23]. Also, supposing 
that kx is the dominant component, a natural dichotomy, 
as observed [25 , would arise in Josephson interferometry, 
between junctions along the face ± x, where the phase 
would be constant across the face, and ± where it 
would alternate on a tight length scale. 

However the energetics that would drive c-axis domain 
formation are not clear: Meissner screening of the kx^iky 
edge currents cancels the long-range fields that encourage 
domain formation in conventional ferromagnets. Also, 
the absence of magnetization at microscopic surface fea- 
tures, such as cleave terraces, implies a tighter limit on 
c-axis domain width than ~400 nm, such that c-axis do- 
mains should be considered a possibility rather than a 
likelihood. 

There is less experimental guidance on the possibil- 
ities for PrOs4Sbi2. Because PrOs4Sbi2 is cubic, if 
the orbital OP turns out to break time reversal sym- 
metry then a more complex domain structure than in 
Sr2Ru04 is expected E.g. in Ref. a singlet OP 
^ {dyz + idxz + Odxy) IS illustratcd, which would give six 
types of domains: ±x, ±y and ±z. If a domain structure 
does exist in PrOs4Sbi2, then the /iSR data indicate that 
the domains would likely be smaller (as scaled by A) than 
in Sr2Ru04 [30]. 

In conclusion, magnetic scans of Sr2Ru04 and 
PrOs4Sbi2 have been presented. No static magnetiza- 
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tion remotely near the TRSB fields observed by /iSR in 
both materials, or theoretical expectation for ± ipy 
order in Sr2Ru04, was observed. 
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